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We present a short overview on the studies of transverse momentum dependent parton distribution functions 
of the nucleon. The aim of such studies is to provide a three dimensional imagining of the nucleon and a 
comprehensive description of semi-inclusive high energy reactions. By comparing with the theoretical frame¬ 
work that we have for the inclusive deep inelastic lepton-nucleon scattering and the one-dimensional imaging 
of the nucleon, we summarize what we need to do in order to construct such a comprehensive theoretical frame¬ 
work for semi-inclusive processes in terms of three dimensional gauge invariant parton distributions. After that, 
we present an overview of what we have already achieved with emphasize on the theoretical framework for 
semi-inclusive reactions in leading order perturbative QCD but with leading and higher twist contributions. We 
summarize in particular the results for the differential cross section and the azimuthal spin asymmetries in terms 
of the gauge invariant transverse momentum dependent parton distribution functions. We also briefly summarize 
the available experimental results on semi-inclusive reactions and parameterizations of transverse momentum 
dependent parton distributions extracted from them and make an outlook for the future studies. 
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I. INTRODUCTION 

With the deeply going of the study of the nucleon structure, 
three dimensional imaging has become the very frontier and 
a hot topic in recent years. It is commonly recognized that 
the three dimensional imaging contains much more abundant 
physics on the nucleon structure and the properties of quantum 
chromodynamics (QCD). The study was initially triggered by 
the experimental finding of striking single-spin asymmetries 
(SSA) in inclusive hadron production in hadron-hadron col¬ 
lisions with transversely polarized hadron [1]. Gradually it 
grows into a field aiming at a comprehensive three dimen¬ 
sional description of the nucleon structure including spin and 
transverse momentum dependences. 

The one dimensional imaging of the nucleon is provided 
by the Parton Distribution Functions (PDFs) such as the num¬ 
ber densities, q(x), the helicity distributions, A q(x), and the 
transversities, 6q(x), for quarks of different flavors in the nu¬ 
cleon. These one dimensional PDFs can be studied in inclu¬ 
sive high energy reactions and are necessary for the descrip¬ 
tion of such inclusive processes. In the three dimensional case, 
i.e. where the parton transverse momentum is also considered, 
not only the direct extensions of these distribution functions 
to include transverse momentum dependences are involved, 
but also many other correlation functions that describe in par¬ 
ticular the correlations between the transverse momenta and 
spins such as the Sivers function, the Boer-Mulders function, 
the pretzelocity etc. exist. They are generally called trans¬ 
verse momentum dependent (TMD) PDFs. Moreover, higher 
twist effects become also important and need to be consid¬ 
ered consistently. The content of the studies is therefore much 
more abundant and more interesting. These TMD PDFs can 
be studied in semi-inclusive reactions and are necessary for 
the description of such processes. 

The study on the three dimensional imaging of the nucleon 
is in a rapid developing phase and it is not so easy to make 
a comprehensive overview of all different aspects of the stud¬ 
ies. Here, we choose to arrange the review in the following 


way: First we will make a brief review of what we did in one 
dimensional case with inclusive deep inelastic lepton-nucleon 
scattering (DIS). In this way, we hope that we can find out the 
main line of what we need to do in three dimensional case. 
After that we will try to summarize the progresses already 
achieved along this line and what we need to do next. Such a 
brief review of the one dimensional case will be presented in 
Sec. 2. In Sec. 3, we will make a short summary of TMDs 
defined via quark-quark correlator. In Sec. 4, we will present 
a brief overview of what we have for constructing the theoreti¬ 
cal framework of semi-inclusive processes. In Sec. 5, we will 
make a short summary of the available experimental results 
and TMD parameterizations extracted from them. Finally we 
will make a short summary of this review in Sec. 6. 

This overview article is an extended version of a plenary 
talk at the 21st international symposium on spin physics 
(Spin2014) [2]. As can be imagined that the simplest and ba¬ 
sic picture is what we have at the leading order in perturba¬ 
tive QCD (pQCD) and at the leading twist. Hence, there are 
also two major directions in theoretical developments towards 
a comprehensive description of the semi-inclusive processes. 
One is to take higher order pQCD into account, and the second 
is to consider higher twist contributions. These contributions 
are important not only for higher accuracy but also for con¬ 
sistency. The major progresses that have been made in recent 
years are also in these two directions separately, i.e. either 
at the leading twist but leading and higher order in pQCD or 
leading order in pQCD but leading and higher twists. The 
talk [2] was mainly concentrated on the second direction. For 
higher order pQCD contributions where evolutions of PDFs 
are involved, an overview talk was also presented by Daniel 
Boer in the same conference [3]. There are also many other 
reviews and monographs (e.g. [4, 6, 7]). The study for higher 
order in pQCD and higher twists seems to be rather difficult 
and even the factorization properties are unclear [5]. In this 
article, we follow the same line as in the talk [2] but briefly 
summarize the progresses in the studies on QCD evolutions 
and refer the interested readers to those reviews. 
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II. INCLUSIVE DIS & THE ONE DIMENSIONAL 
IMAGINING OF THE NUCLEON 


Our studies on the structure of a fast moving nucleon started 
with inclusive DIS such as + N —> e~ + X. We recall 
that, under one photon exchange approximation, the differ¬ 
ential cross section is given by the Lorentz contraction of the 
well-known leptonic tensor IF v (l, A/) and the hadronic ten¬ 
sor Wfiv(q, p, S ), i.e.. 


2a 2 , d 2 l' 

dcr= —^-LX v (l,l',A,)W MV (q,P,S)—. (2.1) 

sQ 4 2 E v 

The leptonic tensor is calculable and is given by, 

LfiviU l', A/) = 2 (lfjl' v + lyl'j — g/jyl ■ /') + HAie^vpa-Pq a ■ (2.2) 

Information on the structure of the nucleon is contained in the 
hadronic tensor defined as, 

Wpviq, P,S) = |^(°)| WIA(0)| P,S) 

Z7T x 

X (2nftf(p + q- px). (2.3) 

Here, l and p denote the 4-momenta of the lepton and the nu¬ 
cleon respectively, those with prime are for the final states; 
A stands for the helicity, and S for the polarization vector 
of the nucleon. We use the light-cone coordinate and define 
the light-cone unit vectors as h — (1,0, Oj.), n - (0,1,0 X ), 
n± = (0,0, ft l), so that a general four-vector can be decom¬ 
posed as = A + W 1 + A~rf + A±, with A ± = (A 0 ± A 3 )/ s/2, 
and A ± = (0,0, A ± ). We work in the center of mass frame of 
the y*N and choose the nucleon’s momentum as z-direction 
so that p and S are decomposed as, 

M 2 

ff = p+rf + — rf, (2.4) 

2 P + 

T) + A/[ 

= A—ff + S^- A - nr (2.5) 

M 7 2p + K J 

The Bjorken variable is defined as xp = Q 1 /2p ■ q, q - 
—; <BP + n + n Q 2 /(2xbP + )', and we also define y - p ■ q/p ■ l. 

The theoretical framework for inclusive DIS has been con¬ 
structed in the following steps. First, we studied the kine¬ 
matics and obtained the general form of the hadronic tensor 
by applying the basic constraints from the general symmetry 
requirements such as Lorentz covariance, gauge invariance, 
parity conservation and Hermiticity, e.g., 


where wf v \q, p) and W^\q, p,S) are given by, 

W$>(q,p) - 2(- iV + ^Vj(x, Q 2 ) 

1 9 

+ —pp(Shi + 2 xp M )(q v + 2xp v )F 2 (x, Q~), (2.10) 

TMs^pcrq 9 
p-q 

x is rr g\(x, Q 2 ) + (S* - —p ,r )g 2 (x, Q 2 )}, (2.11) 

1 p-q 1 

respectively. We found out that the hadronic tensor is deter¬ 
mined by four independent structure functions F\, F 2 , gi and 
g 2 , where the first two describe the unpolarized case and the 
latter two are needed for polarized cases. 

Our knowledge of one dimensional imaging of the nucleon 
starts with the “intuitive parton model” that is very nicely for¬ 
mulated e.g. in [8]. Here, it was argued that, in a fast moving 
frame, because of time dilation, quantum fluctuations such as 
vacuum polarizations can exist quite long. In the infinite mo¬ 
mentum frame, such fluctuations exist forever. In this case, 
a fast moving nucleon can be viewed as a beam of free “par- 
tons”. The probability of the scattering of an electron with a 
nucleon is taken as the incoherent sum of that of the scatter¬ 
ing with each individual parton, more precisely, a convolution 
of the number density of the parton in the nucleon with the 
probability of the scattering with the parton, i.e., 

\M(eN -> eX)\ 2 = ^ f dxf q (x)\M(eq -> eq)\ 2 , (2.12) 

where f q (x) is the number density of parton of flavor q in the 
nucleon. In this way, we obtained the famous results [8], 

F 2 (x, Q 2 ) = 2xFi(x, Q 2 ) = 2 e 2 q xf q (x), (2.13) 

g l (x,Q 2 ) = Yj e l A M x ')’ ( 2 - 14 > 

q 

gi(x, Q 2 ) + g 2 (x, OF) - 2 e\x5f q (x). (2.15) 
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Here, we would like to point out that, with this intuitive parton 
model, we are doing nothing else but the impulse approxima¬ 
tion that we often use in describing a collision process where 
we do the following approximations, 

• during the interaction of the electron with the parton, 
interactions between the partons are neglected; 



q^W qv (q,p,S) = 0, 

(2.6) 

W llY (q,p,-S) = W lv (q,p,S), 

(2.7) 

W* Y (q,p,S) = W w (q,p,S), 

(2.8) 


where A denotes the results of A after space reflection, i.e., 
A 11 - A f ,. The general form of the hadronic tensor is given by 
the sum of a symmetric part and an antisymmetric part, 

W MV (q,p,S) - W%\q,p) + iWfi(q,p,S), (2.9) 


• the electron interacts only with one single parton each 
time; 

• the scatterings of the electron with different partons are 
added incoherently. 

Although the physical picture of the intuitive model is very 
clear and the model is elegant and practical, we are not sat¬ 
isfied with the formulation because it is partly qualitative or 
semi-classical hence it is not easy to control the accuracy. A 
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proper formulation should be based on quantum filed theory 
(QFT) and is obtained by starting with the Feynman diagram 
Fig. 1(a). Here, from this diagram, we obtain immediately 
that, 


W$(q,p,S)= ± J ^Tr[H<°J(k,q)ft°\k,p,S)], 

(2.16) 

where k is the 4-momentum of the parton. 

H$(g, k) = + 4)y v (2n)5 + (.(.k + q) 2 ), (2.17) 

is a calculable hard part. The matrix element 

ft°\k,p,S) = J d 4 ze ikz (p,S\mHz)\p,S), (2.18) 

is known as the quark-quark correlator describing the struc¬ 
ture of the nucleon. By taking the collinear approximation, 
i.e. taking k ~ xp, and neglecting the power suppressed con¬ 
tributions i.e. the o(M/Q) terms, we obtain 


+^) 


q,,q v (q + 2 xp) p (q + 2 xp\ 


q 


2 xp ■ q 


]f q {x). 

(2.19) 


This is exactly the same result as that obtained from Eq. (2.12) 
based on the intuitive parton model. At the same time, we 
obtain the QFT operator expression of f q (x) defined via the 
quark-quark correlator given by Eq. (2.18) as, 

/,(*) = f ^e ixp+z ~(p\m^(z)\p). (2.20) 

By inserting the expanded expression of the field operator if/{z) 
in terms of the plan wave and the creation and/or annihilation 


J 


operators, we see clearly that f q (x) is indeed the number den¬ 
sity of parton in the nucleon. However, from this expression, 
we see also immediately a severe problem, i.e. this expression 
is not (local) gauge invariant! We understand that the physical 
quantity has to be gauge invariant and therefore have to find a 
solution for this. 



FIG. 1: Examples of the Feynman diagram series with multiple gluon 
scattering considered for y* + N —> q + X with (a) j = 0, (b) j = I 
and (c) j = 2 gluons exchanged. 


The gauge invariant formulation is obtained by taking into 
account the multiple gluon scattering shown by the diagram 
series in Fig.l(a-c). This is clear since (local) gauge invari¬ 
ance implies the existence of the gauge interaction that needs 
to be taken into account. In this way, we obtain. 


W/jyiq, p,S) = J] W ( Jy (q, p, S ), (2.21) 

j=0 

where W\j v (q, p,S ) represents the contribution from the dia¬ 
gram with exchange of /-gluon(s). They are all expressed as 
a trace of a calculable hard part and a matrix element depend¬ 
ing on the structure of the nucleon. E.g., corresponding to Fig. 
1(b), we have j = 1, and W p J(q,p,S) is given by. 


Wj»(q,p,S)= J] Wj^\q, p, S ), (2.22) 

c=L,R 

W%'\q,p,S) = ± J ^ ^(ki, k 2 ,q)$\k i, k 2 ,p,S)], (2.23) 

4 >< 'p\k\,k 2 ,p,S) = J ^V' ilZ+te ^ l)y <F,5|0'(O)A p (y) i A(z)| / 7,5), (2.24) 


where c in the superscript represents different cuts (left or right) in the diagram. Similarly, corresponding to Fig. 1(c), we have. 


W%>(g,p,S)= J] W (2 - c \q,p,S), 

c=L,M,R 


W% c \q,p,S) - i- 
(j> ( pl{ki,k 2 ,k,p,S) = 


f 

f 


d 4 k i d 4 k 2 d 4 k 

- r - ~t\ - rTr[Hi; )p (ki ,k^,k,q)d>^(k u k 2 , k, p,S)], 

(2n) 4 (2 tt) 4 (27t) 4 p pc 

d 4 yd 4 y'd 4 z^ y+ik<y - y)+ik ^- y \p, S \m g A p (y)gA a (y')ils{z)\p, S ). 


(2.25) 

(2.26) 
(2.27) 


The matrix element is now a quark-j-gluon(s)-quark correlator. We also immediately see that none of such quark-y-gluon(s)- 
quark correlators is gauge invariant. 

To get the gauge invariant form, we need to apply the collinear expansion proposed in Refs. [9-1 1], which is carried out in the 
following four steps. 
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(1) Make Taylor expansions of all hard parts at kj = x,p, e.g., 

1 d 2 Hfj(x) 


dH {0) (x ) 

mk, q )=m X )+-^^oj n p ‘k 


dk n 


,p kpfta * K ' + ' 


fv p iy/v (T 


H% L)P (k i, k 2 , <?) = - *2) + 


dff p ; L)P (xi,X 2 ) ^ <3f/ pv ’ i)P (xi,X 2 ) ^ 

— (x» kltr' H- T-. - <1) k2o-' + ■ 


dk la 

/ _ 


dk■ 


2 (T 


and so on, where w p p is a projection operator defined by a>p = gp - n p n p . 

(2) Decompose the gluon field into longitudinal and transverse components, i.e., 


A p (y) = A + (y)n p + u p A p ,(y). 


(3) Apply the Ward identities such as, 

- -H«\x x) 
dk p pv (X ’ X) ’ 

(1 ,L) P( 


dH Hr ( X 1 ,X 2 ) _ ~( 2 ,L)p<r, v „ , /V(2 ,M)po-, 


dk] /T 
P P Hpy L)p (x i,x 2 ) 


= -h; 


(xuXux 2 )-mr* K '(xuXux 2 ), 


H$(X t) 


X2 - X\ - l€ 

P P H { p; L)p(T (x i, x, x 2 ) =--- -H^ l)(T (xi , x 2 ). 

X — X\ — 16 

PpHpy M)p,T (x l,X,X 2 )~ 


—H^: L)<T {x i,x 2 )- 


X2 — X\ — 16 

(4) Add all terms with the same hard part together and we obtain, 

W llv (q, p, S ) = ^ W$(q,p,S), 
i 

■A 


X\ — x + le 




W^(q,p,S ) = i- J ±L T r[H^{x) ¥°\k,p,S)], 

d% 


W (P) {q,p,S) 2n 


WV(q,p,S)=± ; f 


(2tt) 4 (2tt) 4 (2tt) 4 j 


^ TrfW p 2 / ;,p,r (xi, x 2 , x)w/ 


(2.28) 

(2.29) 


(2.30) 

(2.31) 

(2.32) 

(2.33) 

(2.34) 

(2.35) 

(2.36) 

(2.37) 

(2.38) 

(2.39) 


where O^’s are the gauge invariant un-integrated quark-quark and quark- /-gluon(s)-quark correlators given by, 

&°\k,p,S) = J d*yeto{p,S\fcO)JXO;yWy)\P.S)' 
a6\k u k 2 ,p,S) = f d*yd*ze^*™(p,S\ftO)W;z)Dp(z)£(z;y)W)lP,S), 

¥£(k u k 2 ,k,P,S) = J AA , ^e ,l4y+, ' % '“ y)+ * 2fe " y \p,5|^(0)£(0;y)D p 0;)i:0;;/)D^0;')i:(/;z)«A(z)|p,5), (2.42) 


(2.40) 

(2.41) 


D(y ) is the covariant derivative defined as D p (y) = -id p + gA p (y). The factor _£(0;y) is obtained during summing different 
contributions with the same hard part together and is given by. 


-£(0;y) = X t (oo;0)X(oo;j’), 

6(oo- y ) - P< r*£orwrto = i - ; g £* 


drA + (y + ,rxx + ( 


/~»oo /'»(X 

~‘ gf Jr Jr l 


(2.43) 


dp A (y ,yj_)A (y ,p ,y ± ) + ■ ■ ■ . 


(2.44) 
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where P stands for the path ordered integral. -£((); y ) is nothing else but the well-known gauge link that makes the quark-quark 
or quark-/-gluon(s)-quarl< correlator, thus also the PDFs defined via them, gauge invariant. 

In this way, we have constructed a theoretical framework for calculating the contributions to the hadronic tensor at the leading 
order (LO) in pQCD but leading as well as higher twists in a systematical way. The results are given in terms of the gauge 
invariant parton distribution and correlation functions (generally referred as PDFs). 


We would like to emphasize in particular the following two further points derived directly from these expressions. 


First, we note that after collinear expansion, the hard parts contained in the expressions for W?’ s such as those given by 
Eqs. (2.37-2.39) are only functions of the longitudinal component x. They are independent of other components of the parton 
momentum k. We can carry out the integration over these components of k’s and simplify them to, 

W$W,S)=-^ 

^ p + dxTi[H?(x)¥°\x,p,S)\, (2.45) 

W«\q,p,S)=± 

fp + dx\p + dx 2 ^ Tr\H ( ? c)p {x \, x 2 )co p <b?(x\,x 2 ,p,S)\, (2.46) 

c=L,R 


p + dx\p + dx 2 p + dx ^ Tr[H^ p(r (xi,X 2 ,x)cL)p a>^ <S>^? (r ,(xi 9 X 2 ,x,p,S)], (2.47) 

c=L,R,M 


where the matrix elements <I>’s are given by, 
d 4 k 


O (0 \x,p,S)= J 2±S(k + -xp + )O m (k,p,S) = J 5’ |iA(0)£(0; y~My-)\p, S ), 

fy 1) (x 1 ,x 2 ,p,S) = J - xip + )5(k2 - x 2 p + )& l) (k u k 2 ,p,S) 

= 5 1-MOmO; z-)D p (z~)£(z-;y-W(y-)\p, S ), 

¥^(x u x 2 ,x,p,S) = J 
J' dy~ dy'~ dz 


(2.48) 


(2.49) 


~ XlP+ ^ k 2 ~ X2 P + )5{k + - xp + )¥ 2 \k u k 2 ,k,p,S) 


2n 2n 2n 


.eWr+W-y-Wninr-r^s \t(OmO;y-)D p (y-)£(y-;ynDAyn£(y'-\Z-)>k(z-)\p,S). (2.50) 


From these expressions, we see explicitly that only x ,-dependences of the quark-quark and/or quark-/-gluon-quark correlators 
are involved. This means that only one dimensional imaging of the nucleon is relevant in inclusive DIS. 

Second, due to the existence of the projection operator co/,' ’s, the hard parts can be further simplified to a great deal. They are 
given by. 


H { py(x) = nW®6(x-x B ), 


H, 


<1X)p (xi,x 2 )w/ = n hpy P u>p 6(xi - x B ), 
2q ■ p 


fiv 


H% L)pa '(Xi,X 2 ,x)u>£ co? = 


H 


MV 
{2,M)po- 


2 n N a)p(r 

An p hT+ " 


—]up oj? 5(x i - x B ), 


(2 q ■ p) 2 pv X2 - xb - is 


(.V|, A'2. X.)OJ. j CO iip V ujp us 

(2 q-pY 


2 7T 


pv 


hpv PCT 0 J p co? 5{x - x B ), 


(2.51) 

(2.52) 

(2.53) 

(2.54) 


where h? - y^y v / p + , h? p - yp$y p i/iy Y , h? pa - p + ypi)y p i/y cr $y Y /2 and Np~ v pcr - q y B y p >/iy CT yv are matrices independent of x,’s. 
We insert them into Eqs. (2.45-2.47) and obtain the simplified expressions for the hadronic tensor as, 


1. 


Wpy (q, p,S) = -Tr [h?^>(x B )\, 


W? L \q,p,S) = ^Tr [h? p co p p '^,\x B )\, 


Aq-p 


W? L \q,p,S) = (2q l p)2 fir[ft? p co? 4p?(x B )\ + Tr[N™ p °co? co? ? 2 ?,(x B )\), 
W <2M> (q,p,S) = 


(2.55) 

(2.56) 

(2.57) 


1 


(2<? ■ p ) 2 


Tr [h^co? co? ?p,?(x B )\, 


(2.58) 
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where, for explicitness, we omit p, S in the arguments of the correlators. These correlators are defined as. 


V ( p\x i) = J 


dx 2 &p\xi,X2,p,S) 


■/ 


p + dy- 


e ixp r (p, S |«0mO; y-)*(y-)\p, S >, 


2 71 

= Jdxidx 2 &™(xux 2 ,x,p,S) = J' ^-^-e ,xp ^ y (p,S\tf/(0)Dp(0)£(0;y~)D ir (y-)i/'(y~)\p,S), 


f dxdx2 . ,x 2 ,x,p,S) 

-^2 -^-l ^ 

= P <**2 p + rfy- pT/g- _ fjp + v - +ffx ,_ j)t>v 

J x 2 - x - it 


■ ie 2n 2n 


(P, S \m-C(0; z-)DJz-)DM~mz-;y-W(y-)\p, S ), 


<Pa L \xi) = J 'dxdx 2 rf<bfl(xux 2 ,x,p,S) = J " , ^-^— e lxp+y (p,S\i//(0)D (0)ZJ o -(0)£(0;y )<A(>' )|p,S). 


(2.59) 

(2.60) 


(2.61) 

(2.62) 


We see explicitly that all the involved components of the 
quark-/-gluon-quark correlators depends only on one single 
parton momentum. This means that only quark-/-gluon-quark 
correlators that depend on one single parton momentum are 
relevant in inclusive DIS. 

We emphasize in particular that the results given by 
Eqs. (2.37-2.39) and their simplified forms given by 
Eqs. (2.55-2.62) including the gauge links are derived in the 
collinear expansion. They are just the sum of the contributions 
from the diagram series shown in Fig. 1. This formalism pro¬ 
vides us a basic theoretical framework for describing inclusive 
DIS at LO pQCD but leading and higher twist contributions 
in terms of gauge invariant PDFs. 

The PDFs are defined in terms of QFT operators via these 
quark-quark correlators by expending them in terms of y- 
matrices and basic Lorentz covariants. For example, for 
<!> (0) (x,p,S), we have. 


at leading twist and have clear probability interpretations, 6 of 
them contribute at twist-3 and the other 3 contribute at twist-4. 
We further note that the three time reversal odd terms e/,(x), 
frix) and li(x) vanish in fact in the one dimensional case. We 
keep them in Eqs. (2.65-2.68) for late comparison with frag¬ 
mentation functions. 

We also see that the PDFs involved here are all scale in¬ 
dependent. This is because we have till now considered only 
the FO pQCD contributions, i.e. the tree diagrams. To go 
to higher order of pQCD, we take the loop diagrams, gluon 
radiations and so on into account. After proper handling of 
these contributions, we obtain the factorized form [6] where 
the PDFs acquire the scale (/-dependence governed by QCD 
evolution equations. In practice, PDFs are parameterized and 
are given in the PDF library (PDFlib). 

In summary, for studying one dimensional imaging of the 
nucleon with inclusive DIS, we take the following steps. 


<i> (0) (x) = i[©f (x) + iy 5 Of s (x) + y“Of (x) 

+ y 5 y a &°\x) + i&* *y sOggtt]. (2.63) 


The basic Lorentz covariants are constructed from p a , n„, S a 
and £ a p P cr. We obtain the following general results, 


Of (x) -Me{x), (2.64) 

Of. (x) =AMe L (x), (2.65) 

M 2 

Of(x) =p + n a fi(x ) + Ms ±ap S P r f T (x ) + — -n a fi(x), (2.66) 

P 

M 2 

3>f (x) =Ap + n a gi L (x) + MS Tagrix) + «o?3l(x), 

P 

(2.67) 


Offx) =p + n [p S Ta]hi T (x) - Ms ±pa h(x) 

M 2 

+ AM n{ p n a \h L (x) + —n^S Ta ]h 3 r(x), 


( 2 . 68 ) 


where s ±pcr = s a p pcr n a nP, and the anti-commutation symbol 
= A p B a -A a B p . The scalar functions/(x)’s, g(x)’s and 
/z(x)’s are the corresponding PDFs. There are totally 12 such 
functions, 3 of them, i.e. /i(x), gir(x) and h\j{x), contribute 


• General symmetry analysis leads to the general form of 
the hadronic tensor and/or the cross section in terms of 
four independent structure functions. 

• Parton model without QCD interaction leads to LO in 
pQCD and leading twist results of structure functions 
in terms of (/-independent PDFs without (local) gauge 
invariance. 

• Parton model with QCD multiple gluon scattering after 
collinear expansion leads to LO in pQCD, leading and 
higher twist contributions in terms of (/-independent 
but gauge invariant PDFs. 

• Parton model with QCD multiple gluon scattering and 
“loop diagram contributions” after collinear approxima¬ 
tion, regularization and renormalization leads to leading 
and higher order pQCD, leading twist contributions in 
factorized forms in terms of (/-evolved and gauge in¬ 
variant PDFs. 

In the following, we will follow these four steps and sum¬ 
marize what we have achieved in the three dimensional case. 
As did in [2], we will mainly focus on the theoretical frame¬ 
work at LO pQCD but taking leading and higher twist contri¬ 
butions into account consistently. Before that, we would like 
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to emphasize the following two of the historical developments 
that may be helpful to us in constructing the theoretical frame¬ 
work for the TMD case. 

First, as mentioned, the study of three dimensional imaging 
of the nucleon was triggered by the experimental observation 
of the single-spin left-right asymmetries (SSA) in the inclu¬ 
sive hadron-hadron collision with transversely polarized pro¬ 
jectile or target. It was known that pQCD leads to negligibly 
small asymmetry for the hard part [ 12] but the observed asym¬ 
metry can be as large as 40% [13]. The hunting for such large 
asymmetries lasts for decades with the following milestones: 

• In 1991, Sivers introduced [14] the asymmetric quark 
distribution in a transversely polarized nucleon that is 
now known as the Sivers function. 

• In 1993, Boros, Liang and Meng proposed [15] a phe¬ 
nomenological model that provides an intuitive physical 
picture showing that the asymmetry arises from the or¬ 
bital angular momenta of quarks and what they called 
“surface effect” caused by the initial or final state inter¬ 
actions. 

• In 1993, Collins published [16] his proof that Sivers 
function has to vanish due to parity and time reversal 
invariance. 

• In 2002, Brodsky, Hwang and Schmidt calculated [17] 
SSA for SIDIS using an explicit example where they 
took the orbital angular momentum of quark and the 
multiple gluon scattering into account. 

• In 2002, immediately after [17], Collins pointed 
out [18] that the multiple gluon scattering is contained 
in the gauge link and that the conclusion of his proof 
in 1993 was incorrect because he forgot the gauge link. 
He further showed that by taking the gauge link into ac¬ 
count the same proof leads to the conclusion that Sivers 
function for DIS and that for Drell-Yan have opposite 
sign. Belitsky, Ji and Yuan resolved [19, 20] the prob¬ 
lem of defining the gauge link for a TMD parton density 


in light-cone gauge where the gauge potential does not 
vanish asymptotically. 

The second historical development that we would like to 
mention concerns the azimuthal asymmetry study in SIDIS. It 
was shown by Georgi and Politzer in 1977 [21] that final state 
gluon radiations lead to azimuthal asymmetries and could be 
used as a “clean test to pQCD”. However, soon after, in 1978, 
it was shown by Cahn [22] that similar asymmetries can also 
be obtained if one includes intrinsic transverse momenta of 
partons. The latter, now named as Cahn effect, though power 
suppressed i.e. at higher twist, can be quite significant and can 
not be neglected since the values of the asymmetries them¬ 
selves are usually not very large. 

The lessons that we learned from these historical develop¬ 
ments are in particular the following two points, i.e., when 
studying TMDs, 

• it is important to take the gauge link into account; 

• higher twist effects can be important. 


Both of them demand that, to describe SIDIS in terms of 
TMDs, we need the proper QFT formulation rather than the 
intuitive parton model. 


III. TMDS DEFINED VIA QUARK-QUARK CORRELATOR 

The TMD PDFs of quarks are defined via the TMD quark- 
quark correlator <I> m, (x, A ± ; p,S) given by Eq. (2.40) (after in¬ 
tegration over k~). A systematical study has been given in 

[23] and a very comprehensive treatment can also be found in 

[24] . Here, we first expand it in terms of y-matrices and ob¬ 
tain a scalar, a pseudo scalar, a vector, an axial-vector and an 
anti-symmetric and space reflection odd tensor part, i.e., 


4> (0) (x, k ± -p, S ) (x, k ± -, p,S) + iysO™ (x, k ± -p,S)+ y a ^\x, k ± - p,S) 

+ y 5 y a ^\x,k ± -p,S) + icx al3 y 5 ®f al} (x,k ± -,p,S)\. (3.1) 

The operator expressions of these coefficients are given by the traces of the quark-quark correlator with the corresponding Dirac 
matrices. For example, for the vector component, we have. 


<$a\x,k ± -,p,S) = iTr[y tt <I> (0) (x,fc ± ;p,S)] = J' dz d 2 z±e 


Kxp+r-kL-tofr s \m-a0; z)fik(z)\p, S ). 


(3.2) 


We then analyze the Lorentz structure of each part by expressing it in terms of possible “basic Lorentz covariants” and scalar 
functions. From <X> (0, (x, Ay; p,S), we obtain the results as [23], 


f A p '» <r 

a>f\x,k ± -p,S) = M\e(x,k ± ) - ±ptT J- T ef(x,k ± )\. 


(x, Aj_; p, S) = M^Ae L (x, k ± ) - k± e T (x, £j_)], 


M 

Ll ' 


(3.3) 

(3.4) 





8 


TABLE I: The 8 leading twist TMD PDFs defined via the quark-quark correlator. A x means that the corresponding term disappears upon 
integrating the quark-quark correlator over cPk ± . 


quark 

polarization 

nucleon 

polarization 

TMD PDFs 

ifX = 1 

integrated over k ± 

name 

U 

U 

T 

fi(x, k±) 
ff T (x, k±) 

0 

fi(x) 

X 

number density 

Sivers function 


L 

gi L (x,k ± ) 


g il(x) 

Helicity distribution 


T 

gf T (x,k±) 


X 

Worm-gear/Trans-helicity distribution 


U 

hf(x, k ± ) 

0 

X 

Boer-Mulders function 

T 

n id 

h 1T (x,k ± ) 


hi T (x) 

transversity distribution 


T(- L) 

hf T (x,k±) 



pretzelosity 


L 

hf L (x,k±) 


X 

Worm-gear/longi-transversity 


®L 0) (*> & l ; P,S) = P + n a [fi(x,k ± ) - ±P<T ^ ST f^ T (x, & l )] + k ±a [f ± (x,k ± ) - — L fr '(*, k ± )] 


e±a P k i: [\Af l L (x, k ± ) - k± ^ T ffr 2 ( x ,k±)\ + ^-n a [f 3 (x,k±) - ±P °^ L T f 3T (x,k±)\, 


M 
k, ■ S 


<b ( a(x, k ± \p,S) — p + n a \Ag\ L (x, k ± ) - -gf T (x, k ± )j + MS Ta g' T (x, k ± ) - s ±a pk ,i _ L g- L (x, k ± ) 


M 

■ k± a [Agi(x,k ± ) - k± gf(x,k ± )\ + ^ -n a [Ag 3L (x,k ± ) - k± ^ T g 3T (x,k±)\, 


^fpaix, k ± \ p,S) - p + ll [p S Ta]hl T (x, k ± ) - 


P + n lp e ±aW k ,i L ± 


M 


hf (x, k ± ) + 


p + n [p k ±a] 

M 


\Ahf L (x, k ± ) - k± lif T (x, fc ± )] 


ST\pk± a ]hf(x, k ± ) - Ms ±pa h(x, k ± ) + Mn [p n lr ]\Ah L (x, k ± ) - k± h T (x, ^j.)] 


M 


M 2 1 „ , ... n^k^ 


1V1 ( 

+ —p[n{pS Ta]h3T(x, k ± ) 


M 


[■ Ahf L (x, k ± ) - k± J’ T hf T (x, A: ± )] ■ 


H\p£±a]pk/^_ 

M 


(3.5) 


(3.6) 


hf(x,k ± )}. (3.7) 


These scalar functions are known as TMD PDFs. There are totally 32 such TMD PDFs. Among them, 8 contribute at leading 
twist and they all have clear probability interpretations such as the number density j\ (x, kj_), the helicity distribution g i / (x, k ± ), 
the transversity Iiit(x, k ±), the Sivers function ff r (x, k ± ), the Boer-Mulders function hf(x, k±) etc.; 16 contribute at twist-3 and 
the other 8 contribute at twist-4. We emphasize that they are all scalar functions of x and k ± , i.e., depending on x and kf. 

If we integrate over d 2 k ± , terms that the basic Lorentz covariants are odd in k ± vanish. Eqs. (3. 3-3. 7) just reduce to the 
corresponding Eqs. (2.64-2.68). At the leading twist, only 3 of 8 survive, i.e. the number density f] (x), the helicity distribution 
gn(x) and the transversity hir(x). 

We show the leading twist TMD PDFs in table I. Those twist-3 TMD PDFs are shown in table II. In these tables, we show 
also the results for the case that X = 1 , i.e. if we neglect the multiple gluon scattering and simply take a nucleon as an ideal 
gas system consisting of quarks and anti-quarks (see e.g. [24]). We also note that the conventions used here have the following 
systematics: /, g, and h are for unpolarized, longitudinally and transversely polarized quarks; the subscript L or T stands for 
longitudinally or transversely polarized nucleon, and those with subscript 1 for leading twist, without number for twist-3 and 
with 3 are for twist-4; the _L in the superscript denotes that the corresponding basic Lorentz covariant is k ± dependent. 

Higher twist TMD PDFs are also defined via quark-/-gluon(s)-quark correlators such as those given by Eqs. (2.59-2.62). 
Many of them are, however, not independent since they are related to those defined via the quark-quark correlator through the 
QCD equation of motion y ■ D(z)i//(z ) = 0. We can get the relations such as, 


(*, k ± ;p,S) = -^ |Pe^(x, k ± ; p, S) + lmip^l(x, k ± ; p, S)], 


x(F™'(x, k ± ;p,S) = -— [Re^(x, k ± ; p, S) + e^Im<p ( Jl(x, k ± ;p,S)]. 


(3.8) 

(3.9) 


It is interesting to see that [35], although not generally proved, all the twist-3 TMD PDFs that are defined via quark-gluon-quark 
correlator ipp ] and involved in SIDIS are replaced by those defined via quark-quark correlator <I> m k 
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TABLE II: The 16 twist-3 TMD PDFs defined via the quark-quark correlator. A x means that the corresponding term disappears upon 
integrating the quark-quark correlator over cPk ± . 


quark 

polarization 

nucleon 

polarization 

TMD PDFs 


if£= 1 

integrated over k ± 

U 

U 

e(x, Cl), 

/■*■(*, Cl) 


0, fi(x,k ± )/x 

e(x), x 

T 

e P (x,k ± ), 

ff l (x,k ± ), 

f P 2 (x,k ± ) 

O 

O 

O 

XXX 


L 

e L (x,k ± ), 

g£(x,k ± ) 


0, g l (x,k ± )/x 

X, X 

L 

T 

e T (x,k ± ), 

g F (x,k ± ) 

0, 0, gir(x, Cl)/* 

x gr(x) 

g’ T (x,k ± ), 


U 

h(x. Cl) 



0 

X 

T 

W (II) 

h p { X, Cl) 



hf T (x,k ± )/x 

X 


T{±) 

h F (x, Cl) 



h lT (x, k ± )/x + k 2 L hf T {x, k ± )/M 2 x 

X 


L 

h L (x,k ± ) 



k\h\ L (x, k ± )/M 2 x 

h L {x) 

U 

L 

/l CtCl) 



0 

X 

L 

U 

g ± (x. Cl) 



0 

X 


We would like to emphasize that fragmentation is just conjugate to parton distribution. A systematic study for the gen¬ 
eral structure of fragmentation function (FF) defined via the corresponding quark-quark correlator is presented in [26]. We 
should have one to one correspondence between TMD PDFs and TMD FFs. E.g., corresponding to the quark-quark correlator 
<f> (0, (k, p, S ) given by Eq. (2.40) and the expanded form Eq. (3.1), we have. 


* {0) (k F ,p,S) 


■sZf 


0|X + (0, oo)tmhX)(hX\fa)£(t, oo)|0>. 


(3.10) 


F <0) (z, k F± -, p,S ) =^[H^ 0) (z, k ± \p,S ) + iy 5 E^(z, k ± \p,S ) + y“E^ 0) (z, k ± \p,S ) 

+ 757 a ^a ) (z,k ± -,p,S) + k ± \p,S )]. 


(3.11) 


For spin-1/2 hadron, we have perfect one to one correspondence to those given by Eqs. (3. 3-3. 7) for parton distributions in the 
nucleon, i.e.. 


zZf(z,k F± -,p,S) = M[E(z,k F± )- 


s± P o-k p F± S% 


Ej(z,k F x )], 


M 

zZpl(z,k F j_;p,S) = M[AE L (z, k F± ) + kp± S T E T (z, k F± )\, 


zH^ 0) (z, k F± -, p,S) = p + n a [D l (z, k F± ) 

T £±apk?pj_ (z, k F j_) 4 ^ 

zE^ 0) (z, k F _ L ; p, S) = p + n a \AG\ L (z, k F± ) + — 


M 

S±pcrk P F± Sj 

h M 
k F ± ■ S t 


Df r (z, kfj_)J + k F±a D ± (z, k F _ l) + Me ±ap S P T D T (z, k F± ), 

i i M~ r £±po'kf F i S 7. t 

D$(z,k F± )\ + —n a [D 3 (z,k F± ) + °UZ’k F± )], 

-Gf T (z, kfj.)] + MS Ta G T (z, k F± ) + E ±a/ 3k/ 3 F _ L G ± (z, k FA _) 


k F ± a \^G F (z,k F j_) ■ 


M 

k F \ ■ St . , . 1 M~ 1 . „ , . . k F± ■ S i 


M 


-G F (z,k F j.)j + — -n a ^AG 3 L(z, k F± ) 


zSy] a (z, k F± ; p,S) - p + n [p S Ta ]HiT(z, k F± ) 

+ 


qk^t, 


-Hf(z,k F± ) + 


M 

p + n[pk F±a] 


■G 3 t(z, kfj_)J, 


[AHf L (z, k F± ) ■ 


Tpa ^^ r , , r - M M L M 

S T[pk F ±a]H F (z, k Fx ) + Ms ±po H(z, k F± ) + [ MAH L (z , k F± ) + • 5 T H' F (z, k F± )j 

kH 


(3.12) 

(3.13) 

(3.14) 

(3.15) 

-//f r (Z, Cfj.)] 


M 2 , 

+ —p\ n [pS Ta]H 3F (z, k F± ) 


M 


-H 3 (z, k F± ) ■ 


ti[pk F ±a] r_ „ 1 . . , k F ± • S 


M 


■| AHf L (z,k F± )- 


FJ- ' 0 T „_L 


M 


^ r (z,k F± )\). (3.16) 


Comparing them with the results given by Eqs. (3. 3-3. 7), we see clearly the one to one correspondence between FFs and PDFs. 
As an example, we show the 8 leading twist components in table III. We do not show the results for the case of X = 1 for FFs. 
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TABLE III: The 8 leading twist TMD FFs for spin-1/2 hadrons defined via the quark-quark correlator. A x means that the corresponding term 
disappears upon integrating the quark-quark correlator over tfk F± . 


quark 

polarization 

hadron 

polarization 

TMD FFs 

integrated over k F± 

name 


U 

Di(z,k F± ) 

Di(z) 


number density 

U 

T 

Df T (z,k F± ) 

X 




L 

G lL (z,k F± ) 

G lL (z ) 


spin transfer (longitudinal) 

L 

T 

Gf T {z,k F± ) 

X 




U 

Hf(z, k F± ) 

X 


Collins function 

T 

7X11) 

Hir(z,k F± ) 

Hit(z) 


spin transfer (transverse) 


n±) 

Hf T (z,k F± ) 




L 

Hf L (z,k F± ) 

X 




This is because even if we neglect the multiple gluon scattering that leads to the gauge link, final state interactions can still exist 
between h and X. In this case, time reversal invariance does not lead to zero results for the T-odd amplitudes. 

For spin-1 hadrons, the polarization is described by the polarization vector S and also the polarization tensor T (see e.g. 
[25] and [26]). The tensor polarization part has five independent components. They are given by a Lorentz scalar Sll , a 
Lorentz vector S^ T = (0, S X LT , S y LT , 0) and a Lorentz tensor S^ T that has two independent non-zero components S XX T and Sj T 
in the rest frame of the hadron. These polarization parameters can be related to the probabilities for the particles in different 
spin states [25]. In this case, the TMD quark-quark correlator E ,0) (z, k F± \p,S ) is decomposed into a spin independent part, 
a vector polarization dependent part and a tensor polarization dependent part, i.e. E (0) (z, k F± \p,S ) = E ,/(0) (z, k Fi _\p,S) + 
E y(0) (z, k F± \ p,S) + E r(0) (z, k F± \p,S ). The spin independent and vector polarization dependent part E c/+y(0) (z, k F± \p,S ) takes 
exactly the same decomposition as that for spin-1/2 hadron given by Eqs. (3.12-3.16). The tensor polarization dependent part is 
presented in [26] and is given by. 


zS£ ( 0 ) (z, k F± ; p,S) = m[s ll E ll (z, k F± ) 


M 


-Et T (z, k F _ l) ■ 


M 2 


-Ej T (z,k F _ l)], 




JcfSl 


k F ±',p,S) - M[^—E' L ^z,k F± ) + - —E' rr (z.k F . )|. 


k F± -p,S) = p + n cr [s LL Di LL (z,k F± ) + — - - ^ LT Df LT (z, k FA _) + kf 


M 


M 

S LT _ ± , , , k 

D lt (z, k F± ) 


S F t • k F± 

~W- 


n - 1 
u \TT 


(z,k F± )\ 


+ kF±a llD ll (z, k F _l) ■ 

+ MSltoDlt(z, k F j _) + TT pa E>' FT {z, k F ±) 


F± ' OTT ' K-F± 


M 2 
k F 1 • S 7 


Dj T {z, k F± ) | 


M r r. / 7 , k Fl _ ■ S LT „_L , , . k F± ■ S TT ■ k Fl _ , I 

+ — ~na[S llE>ill(z, k F± ) H- — - E> iLT (z,k F± ) H-- E> iTT (z,k F± ) J, 


zSl (P \z,k F± ;p,S) = p + n, 




kpj_S l 


: p£T 


M 


p-,_L , , \ . s ±k F± pkF±a-S TT 

-tr lLT (Z,KF±) + m 2 GITT 


(z> k F ±)\ 


+ Sj.pi 


r F± \s uGfj^kpj 


VfJ. 


M 


~G it (z, k F± ) + 

p°- n'-L , 


V F± 


M 2 


-G FT (z,k F _ L )] 


+ Ms±p a S F jG FF {zi^ k F i ) + £±apk F i_Q-S FF G k F ±) 


M~ r 
+ —n a y 


G VJ (z, k F j_) ■ 


s F k F1 pk 


z" n0) (z 


k F ±',p,S ) 


pc l m 

P kl[p£±a](rkpj_ 


FJ.crO TT , 

-Lr, 


M 


\sLLHf LL (z, k F± ) 


M 1 

k F ± 


(z,k F± )], 


M 


Slt TT _ L , , , k F 

EI IL t(z, k F± ) -(- 


(3.17) 

(3.18) 


(3.19) 


(3.20) 


M 2 


-H\tt(z, &f±)J 


+ p + n[ p s ±a ] (r S < [ T Hiijr(z, k F± ) + 

k 


Ms ±pa yS llH L l(z, k F± ) ■ 


P fl[p^±a]crk F ±sS FF 

M 

Sfi.,, , , , k, 

H ft (z, k F _ l) + 


| y / (Z, k F j_) 

S TT ' kp 


M 


M 2 


-Hj r (z„ k F± )] 
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TABLE IV: The 10 tensor polarization dependent TMD FFs for spin-lhadrons defined via the quark-quark correlator. A x means that the 
corresponding term disappears upon integrating the quark-quark correlator over cPk F± . 


quark 

polarization 

hadron 

polarization 

TMD FFs 

integrated over k F± name 


LL 

k F± ) 

D lLL (z ) spin alignment 

U 

LT 

D IL t( Z ’ ^ F± ) 

X 


TT 

D lTT (z, 

X 


LT 

G i/j-fo ^F±) 

X 

L 

TT 

X 

G lTT (z, 


LL 

Hf LL (z,k F G 

X 

T 

LT 

HiL T (z,k F± ), Hf LT (z,k F± ) 

Hilt(z) 


TT 

H lTT (z, kfj_\ H ^f_l) 

X 


T k F ±S L TTi>-L/ l \ , E -kk F± a-kFl.sSj 6 T 
|e_L H LT (z,k F± ) + - 


M 




M~ i n[pE±a]o-kp ± r , k F ± ~ Slt Tr ± . , . k F± ■ Stt • k F ± 

+ ~r{ -- [s LLH 3LL (z,k F± ) +-—- H 3lt (z, k F± ) + 


M 2 


[ 37T 


(Z,&F.l)] 


+ n[pS±a]a-S‘[ T H3iT(z, k F± ) + 


tfi /-j -j 


M 


H 


TiTT( Z ’k F ±)}- 


(3.21) 


We see that, for the vector polarization dependent part, similar to nucleon TMD PDFs, there are totally 32 components, 8 
contributes at leading twist, 16 at twist-3 and the other 8 at twist-4. For the tensor polarization dependent part, there are totally 
40 components, where 10 contribute at leading twist, 20 at twist-3 and the other 10 at twist-4. In Table IV, we list the twist-2 
components for the tensor polarization dependent part. 

If we integrate over d 2 k F± , we have, corresponding to Eqs. (3.12-3.16), for the spin independent and vector polarization 
dependent part. 


-U+V(0) 

,~t/+n o) 


(z',p,S) = ME(z), 
(z;p,S) = AME l (z ), 


, 3 t/+v ( 0 )(, ;/? ,£) _ p+ na o t (z) + Ms ±ap S p D T (z) + —~n a Di(z), 

P + 

M 2 

zE^ +v(0) (z; p,S ) - Ap + n a G 1L (z ) + MS Ta G T (z ) + A—n a G 3L (z), 


zS^J (0) (z; p, S) = p + n [p STa\H\ T {z) + Ms ±pa H(z) + AMn\j,n„\H L (z) + — -n^Sr a:\HYriz). 


while for the tensor polarization dependent part, we have, 

zE T s ( -°\z-,p,S) = MS LL E LL (z), 
za T p f\z;p,S) = 0, 

M 2 

zE-l < 0 > (z\ p, S ) = p + n a S LL D\ LL (z) + MS LTa D LT (z,k F± ) + —n a S LL D iLL (z), 

P 

zZ T a ( 0 \z-,p,S) = Ms ±pa S p LT G LT (z), 

—T(n\ _ M 2 

Z , — t 7 pa (.ZiP,S') — p [j- H \ u'(z) + Ms ±pa S llMukz) + ~Zl[p£±a]o'Suptl'iLT^Z)' 


(3.22) 

(3.23) 

(3.24) 

(3.25) 

(3.26) 


(3.27) 

(3.28) 

(3.29) 

(3.30) 

(3.31) 


We see that, for the spin independent and vector polarization dependent parts, 12 components survive, 3 of them contribute 
at twist-2, 6 at twist-3 and the other 3 at twist-4. This is exactly the same as those for PDFs for nucleon and we have exact 
one to one correspondence between the results given by Eqs. (3.22-3.26) and those given by Eqs. (2.64-2.68). For the tensor 
polarization dependent part, there are only 8 components survive, 2 of them contribute at twist-2, 4 at twist-3 and the other 2 at 
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twist-4. This corresponds to the situation for PDFs for vector mesons. We should have a one to one correspondence between the 
tensor polarization dependent FFs for production of spin-1 hadron to those PDFs for spin-1 hadrons. We also listed the twist-2 
components in table IV. 


IV. ACCESSING THE TMDS IN HIGH ENERGY REACTIONS 


The TMDs can be studied in semi-inclusive high energy reactions such as SIDIS e ' + A' —> e +h +X, semi-inclusive Drell-Yan 
h + h —> l + + l~ + X, and semi-inclusive hadron production in <? + e~-annihilation e 4 + e~ —> /? | + lu + X. With SIDIS, we study 
TMD PDFs and TMD FFs, while with Drell-Yan and e + e annihilation, we study TMD PDFs and TMD FFs separately. We now 
follow the same steps as those for inclusive DIS and briefly summarize what we already have in constructing the corresponding 
theoretical framework. 

(I) The general forms of hadronic tensors: For all three classes of processes, the general forms of hadronic tensors have been 
studied and obtained. For SIDIS, it has been discussed in [27-30] and it has been shown that one need 18 independent structure 
functions for spinless h. For Drell-Yan, a comprehensive study was made in [31] and the number of independent structure 
functions is 48 for hadrons with spin 1/2. For e + e~-annihilation, the study was presented in [32] and one needs 72 for spin-1/2 
h\ and /zo - The results are systematically presented in these papers and we will not repeat them here. However, we would like to 
present as an example for the general form of the differential cross section for e N —> e~hX. It is given by, 


dcr 


a;, 


dxdydzdi//d 2 pi 1± 

„2 


xyQ- 


(l + ~)(V//(/ + A/Tlu + ATul + AiAT'll + S ±Tut + AiS jTlt), 


Tuu - y —— ( Fuu,t + sFuu,l + V 2e 0 +e)F c °^ 
v 2 

Tul = y 77 ^( V 2s (! + s)F S y L *’' sin 0,, + sF s ^ 4h sin 20,,), 


rCOS 07 , / j-iCOS 20 /, \ 

vu COS 0;, + sF uv ” COS 20,, j, 


v 2 

Tlu = — V 2s (l - e)F*£* h sin0,„ 

2 

Tll = —{ Vl - £2p LL + V2s(l - e)F c L °^ cos 0 A ), 

T,jt = V 2£ ( ] + «)F,7/ S sin0 s + + eF^^ s) ) sin(0 A - 0 S ) 

+ sF'™ (0 ' ,+0s, sin(0,, + 0s) + y/2s{ 1 + s)F s J'j 2<f ’ l ^ 4 ’ s) sin(20 A - <p s ) + sin(30;, - 0s)], 

Tlt = y2_[V2e(l-e)F“ s0s cos 0 S + Vl - cos(0 A - 0s) + ^2e(.\ - e)F™ (2H ~* s) cos(2$ h - 0 S )], 


(4.1) 

(4.2) 

(4.3) 

(4.4) 

(4.5) 

(4.6) 

(4.7) 


where s = (1 - y - \y 2 y 2 )/( 1 - y + ^y 2 + ^y 2 y 2 ), y = 2Mx/Q\ the azimuthal angle 0 is that of the out going lepton /' around 
the incident lepton beam with respect to an arbitrary fixed direction, which in case of transversely polarized target is taken as the 
direction of In the deep inelastic limit, neglecting power suppressed terms, dtjj — d(ps- 

From Eqs. (4. 1 -4.7), we see explicitly that the 18 structure functions F’s are determined by the different azimuthal asymmetries 
in different polarized cases. These different azimuthal asymmetries are just defined by the average value of the corresponding 
trigonometric functions. E.g., 


* s) = (sin(0 A - <t>s))ur 
A f? h+ * s) = <sin(0„ + <Ps))ut 


psin pSin 

r UT,T + UT,L 

2(Fuu,t + eFuu,l) 

psin (<Ph+<Ps) 

UT _ 

2{F uuj + sF UUtL ) 


(4.8) 

(4.9) 


We also like to emphasize that they are the general forms independent of parton model and are valid at leading and higher 
twist and also leading and higher order in pQCD. 

(II) LO in pQCD and leading twist parton model results: These are the simplest parton model results and can be obtained 
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easily. E.g., for SIDIS, 
daW 


-(t' 

1 V 1 


( 0 ) , 1 qH 0 ) 
UV + ' 


dxdydzd(f>sd 2 pi, 1 _ xyQ 2 

+ 2i 

- y)C[wihf L Hf] sin(2 <p h ). 


ffl = A(y)C[f\D\] + 2(1 - y)C[ Wl hfHf] cos(2/,,), 


LU 
J-u-L 


A K°i 


■ A,AT^ +S ± T}Z + 


r(0) 
' J- 7 UT 


A,s ± r%>), 


K°l = 2 d 

T?u = °> 


= C(y)C[gi L £»i], 

Tyj = A(y)C[w 2 ff~ T Di ] sin(0/, - / s ) + 2(1 - yX^v^/iiT-// 1 ] sin(/ ; , + / s ) + 2(1 - >')C[vv 4 /?^// | L ] sin(3//, - (p s ), 
= C(>')C[-w 2 gr/ /J|] cos(/ ; , - / s ), 

where A(y) = 1 + (1 - y) 2 , C(y) = y(2 - y), and C[vv, / /)| denotes the convolution of / and D weighted by w,-, i.e., 


C[wifD] 


= x ^ e 2 J' d 2 k ± d 2 kf ± d (2) (k 


±-k Fj _- Ph T /z)wi(k ± , k F± ,p hT )f q (x, k±)D q ~ >hx (z, k Fx ), 

I- 


(4.10) 

(4.11) 

(4.12) 

(4.13) 

(4.14) 

(4.15) 

(4.16) 


(4.17) 


where the weights w, ’s are given by, 


wi(k x ,k F± ) = 


-2(Pi,t ■ k F± )(phT ■ k ± ) + (k ± ■ k F± ) 


MM h 


w 2 (k ± ,k F± ) = - 
w 2 (k ± ,k F± ) = - 


PhT ■ k x 
M ’ 

PhT • k F j_ 
M h ’ 


{phT ■ ^j_)(^j. ■ k F j_) + k-^iphT ■ k F j_) 

M 2 M h 

_ 2(// i7 - ■ k ± ) 2 (pi, T ■ k F j_) 

M 2 M h 


w 4 (k ± ,k F± ) = 


(4.18) 

(4.19) 

(4.20) 


(4.21) 


where p F r = Phrl\PhT\ is the corresponding unit vector. The 
result can be obtained from those given e.g. in [30] by ne¬ 
glecting all the power suppressed contributions. 

From Eqs. (4.10-4.16), we see in particular that, at leading 
twist, there exist 6 non-zero azimuthal asymmetries in differ¬ 
ent polarized cases, i.e., 


<cos2«® = U 

V uu A(y) C[f\Dx] 

(4.22) 

a , x( 0) (1 -y)C[wihf L Hf] 

UL A(y) C[f\D\] 

(4.23) 

, WO) C[w 2 ff- T D X ] 

<sm(/,, <Ps)) UT - 2C[fiDi] ’ 

(4.24) 

VS ” UT A(y) C[f\D\] 

(4.25) 

x w(0) ( l-y)C[w 4 hf T Hf] 

<sm(30 /7 cf> s )) UT A(y) C[flDl] > 

(4.26) 

1 


<COS(/;, - (f> S ))fl = 


C(y) C[-w 2 g\rP\] 
2A(y) ClfiDy] 


(4.27) 


and they are determined by Boer-Mulders function hf con¬ 
voluted with Collins function Hf, the Worm-gear (longi- 
transversity) hf L convoluted with Collins function Hf, the 
Sivers function ff r convoluted with l )\, the transversity h \ T 
convoluted with Collins function Hf, the Worm-gear (trans- 
helicity distribution) gf T convoluted with Collins function 
Hf. The azimuthal asymmetries Afff^’ 1 ’^ are due to Sivers 
and Collins effects and are often referred as Sivers asymmetry 
and Collins asymmetry respectively. 

We would like to emphasize that the results given by 
Eqs. (4.10-4.21) is a complete parton model result at LO in 
pQCD and leading twist. It can be used to extract the TMDs 
at this order. Any attempt to go beyond LO in pQCD or to 
consider higher twists needs to go beyond this expression. 

(Ill) LO in pQCD, leading and higher twist results: For the 
semi-inclusive processes where only one hadron is involved, 
either in the initial or the final state, it has been shown[33-37] 
that the collinear expansion can be applied. Such processes in¬ 
clude: semi-inclusive DIS e~+N —> e~ +q{ jet)+ X, ande + e~- 
annihilation e' +e —> h+q(jet)+X. By applying the collinear 
expansion, we have constructed the theoretical frameworks 
for these processes with which leading as well as higher twist 
contributions can be calculated in a systematical way to LO in 
pQCD. The complete results up to twist-3 have been obtained 
in Refs. [35-37]. For polarized e~ + N —> e~ + q(jet) + X, the 
simplified expressions for the hadronic tensor are very similar 
to those for the inclusive DIS given by Eqs. (2.55-2.58), 


Wfff si \q,p,S,k x ) = ^Tr 


(4.28) 


















14 


W$ L ’ si \q,p,S,k ± ) 


1 


Tr[/? pv )p w p p (Pp; L \x B ,k ± )], 


4 q ■ P 
1 


Wfi L ’ s %, p,S,k ± ) = _ 2 {Tr[h ( ^ p (v p p ^, L \x B , kj] + Tr[lV P y )ptr w p p a>f <pf£(x B , k±)]}, 


W% M ’ s 'Xq,p,S,ki_) = 


(2 q ■ p) 2 

1 


(2q ■ p) 2 *" L pv p 

and the complete results up to twist-3 are given by, 

„2 „2 


Tr[/i p y' pcr w p p ojf kj]. 


dcr 


2na z em e q 


dxdyd 2 k ± Q 2 y 


XWuu + AfW lu + S j fWuT + AAVul + AiA'Wll + A/S ±4 Vlt), 


2x\k I 

'Wuuix, k ± , (p) = A(y)fq(x, k_ L ) - B (y)f^(x , k ± ) cos cf>. 


2x\k I 

AV LU (x,k ± ,(f>) = - ^—^—D(y)g ± (x, k ± ) sin </>, 

nV(/ r (^,0, </> s ) = ^A(y)/ 1P (x, k ± ) sin(0 - 0 S ) 

2xM k 2 

+ *J sin(20 -0s)- fr(x, k ± ) sin 0 S }, 

*Wul(x, kj_,</>) - - 2X Q^ B(y)f£ (x, k±) sin 0, 

2 x \k I 

nV LL (x, <f>) = C(y)g\i(x, k ± ) - ~ ^ £>(y)g p (x, kj cos 

AV Lr (x, k ± , 0, 0 S ) = ^C(y)gfj.(x, * x ) cos(0 - 0 S ) 

2jcM k^ 

~ —^~D(y)[g T (x, kj cos 0s - ^4_g p (x,k_L)cos(20 - 0 S )]. 


(4.29) 

(4.30) 

(4.31) 

(4.32) 

(4.33) 

(4.34) 

(4.35) 

(4.36) 

(4.37) 

(4.38) 


where B(y) = 2(2 - y) y/1 - y, D(y) = 2y yj 1 - y. For unpolarized e + N —> e + q(jet) + X, the results up to twist-4 have also 
been obtained [34], 


dcr uu 
dxdyd 2 k± 


<r em e 2 q 

Q 2 y 


\My)f\{x, kj - 2 B(y)^xf ± (x, k ± ) cos 0 - 4(1 - y)^-x[^ lu (x, k ± ) - 
+ 8(1 - y)—^~Mx, k ± ) - 2A(y)^-x[(/3 t 3 2 ’ i, - L (x, k ± ) - <pf’ L)± (x, kj]}. 


0 3 )_L (x, kj_)]cos 20 

(4.39) 


These results are expressed in terms of the gauge invariant TMD PDFs or FFs and can be used as the basis for measuring these 
TMDs via the corresponding process at the LO in pQCD. 

We would like in particular to draw the attention to the results for e + + e~ —» h + q(jet)+X for h with different spins [37]. Here, 
for hadronic tensor, we obtain again very much similar formulae also for this process, e.g., corresponding to Eqs. (4.28-4.30), 
we have. 


Wf’ sl Xq,p,S,k’ ± \e + e-) = ^Tr [h^Z {0 \z B ,k' ± )\, 


Wit”\q, p,S,k' ± \e + e ) = Tr [h^ p (Op e£! \z B , k' ± )\, 

t+p • q 

1 


W^'Xq, p,S,k’ ± \e + e-) = ^Tr[^w/ k'J + N$ p *cjf cofS™(z B , k ’ ± )], 

W^ v Msl \q,p,S,k' ± \e + e~) = -^-^-^Tr *™(z B ,k? ± )\. 


(4.40) 

(4.41) 

(4.42) 

(4.43) 


A complete twist-3 results for differential cross sections, az- for hadrons with spin-0, 1/2 and 1 in [37]. We see in partic- 
imuthal asymmetries, and polarizations have been obtained 
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ular for spin-1 hadrons, tensor polarization is involved, even 
at the leading twist level, we have, for e + e annihilation at the 
Z°-pole, 


ZqT^{y)DiLL(z,p T ) 
2 'Lq T o(y) D l(Z’P T )’ 
y 

LT (y,z,p T ) = 


S?l(y,z,p T ) 


s >m 


sfj{y,z,p T ) 


S z, p T ) 


S n jj ] (y,z,pT) 


2\pr\ 2, 

, P q (y)To(y)Gf LT (z, p T ) 

3 zM 

E ? T q (y)D\(z, Pt) 

2\p T \ 2, 

T o(y) D f L A z ’ pt) 

3 zM 2 

q T q ) {y)D\{z, P T) ’ 

2\p T \ 2 E 9 T q (y)Df TT (z, p T ) 

3 M 2 

Tig To(y)Di(y,p T ) 

2\p T \ 2 2, 

P q (y)Tl(y)G\ TT (z, p T ) 

3 M 2 

2t/ Po(y)Di(y> Pt) 


(4.44) 

(4.45) 

(4.46) 

(4.47) 

(4.48) 


where n and t denote the two transverse directions of the pro¬ 
duced vector meson, one is normal to and the other is inside to 
the production plane. The coefficient Tfj(y) = c[c\[{\ - y) 2 + 
y 2 ] - c\c\[ 1 - 2 y], c* = (c e y) 2 + (c«) 2 and ^ = 2c^; and y 
in this reaction is defined as y = l\lk + . P q (y) = Ti(y)/To(y) 
is the polarization of the quark produced at the Z°-decay and 
T q (y) = —c^Cj[(l -y) 2 +y 2 ] + CjCj[1 - 2y], This is a situation 
that is much less explored till now and is worthwhile for many 
further studies. 

For the above-mentioned three kinds of semi-inclusive pro¬ 
cesses, there are always two hadrons involved. Collinear ex¬ 
pansion has not been proved how to apply for such processes. 
It is unclear how one can calculate leading and higher twist 
contributions in a systematical way. Nevertheless, twist-3 cal¬ 
culations that have been carried out for these processes [38- 
41], practically in the following steps: 

(i) draw Feynman diagrams with multiple gluon scattering 
to the order of one gluon exchange, 

(ii) insert the gauge link in the correlator wherever needed 
to make it gauge invariant, 

(iii) carry out calculations to the order 1 /Q. 

Although not proved, it is interesting to see that the results 
obtained this way reduce exactly to those obtained in the cor¬ 
responding simplified cases where collinear expansion is ap¬ 
plied if we take the corresponding fragmentation functions as 
h-fu notions. 

(IV) TMD factorization and evolution: To describe the 
semi-inclusive high energy reactions mentioned above in 
terms of QCD and parton model, TMDs are needed and the 
factorization theorem has to involve transverse momentum de¬ 
pendence. TMD factorization theorem has been established at 
the leading twist for semi-inclusive processes [42-49]. TMD 
evolution theory is also developing very fast [50-63]. There 
was a dedicated overviews by Daniel Boer [3] in Spin2014. 
There is a dedicated annual workshop series since 2012. We 
refer the interested readers to these talks and overviews. 


V. AVAILABLE DATA AND PARAMETERIZATIONS 

Experiments have been carried out for all three kinds of 
semi-inclusive reactions. The results are summarized e.g. in 
a number of plenary talks in Spin2014 by Marcin Stolarski 
and Armine Rostomyan [64, 65]. Here, we will just briefly 
summarize the main data available and then try to sort out the 
TMD parameterizations that we already have. 

For SIDIS, there are measurements carried out by HER¬ 
MES Collaboration [66-70] at DESY, COMPASS Collabora¬ 
tion [71-78] at CERN, CLAS [79, 80] and Hall A Collabo¬ 
ration [81-84] at Jefferson Laboratory. We list these SIDIS 
experiments in table V and briefly summarize the results in 
the following. 

At DESY, the first measurement on single-spin asymme¬ 
tries for SIDIS with longitudinally polarized target was car¬ 
ried out by HERMES [66] for production of charged pions; 
then for the first time with transversely polarized target in 
[67]. They found non zero Sivers and Collins asymmetries 
(sin(0;, - <Ps))ut and (sin(0/, + (f>s))uT■ Measurements have 
then also carried out for^° and Kaons [68, 69] and also for az¬ 
imuthal asymmetries (cos (f>i,)uu and (cos(20/,))j/£/ in the un¬ 
polarized case [70], 

At CERN, COMPASS has carried out measurements on the 
Sivers and Collins asymmetries in reactions with Deuteron or 
proton targets for production of charged hadrons, pions and 
Kaons [71-78], and also (cos <ph)uu and (cos(2 (j>h))uu in the 
unpolarized case [78]. 

At JLab, CLAS has carried out the measurements [79, 80] 
on (sin(20/,)){/£ for pions with different charges and (sin <ph)w 
for 7T°. Hall A Collaboration has made the measurements [8 1- 
84] on Collins and Sivers asymmetries for tP and K ± , 
(cos (cpi, - 4> s ))lt for iP and (sin(30/, - (f> s ))uT- They are all 
summarized in table V. 

Besides the data from SIDIS, we have now also measure¬ 
ments on the azimuthal asymmetries in e + e~ —> n + n + X 
by Belle [85-87] and Babar collaboration [88], and also pre¬ 
liminary results from BES [89]. For Drell-Yan, there are data 
available on azimuthal asymmetries in e.g. reactions using 
pion beam [90-93], and pp or pD collisions [94, 95]. 

Although the data are still far from abundant enough to give 
a precise control of the TMDs involved, there are already dif¬ 
ferent sets of TMD parameterizations extracted from them. 
We briefly sort them out in the following. 

The first part concerns what people called “the first phase 
parameterizations”, i.e. TMD parameterizations without QCD 
evolutions. Here, we have in particular the following results 
available. We emphasize once more that all the results includ¬ 
ing the figures are taken from these references [96-112]. The 
interested readers are referred to these references for more de¬ 
tails. 

(1) Transverse momentum dependence: This is usually 
taken as [96-100] a Gaussian in a factorized form indepen¬ 
dent of the longitudinal variable z or x, e.g., 

fi(x,k ± ) =/i(x)e _ ^ /< ^ > /7r(^), (5.1) 

Di(z,k F± ) = Di(z)e- t M<&J/„<j? F± )' (5.2) 
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TABLE V: Available measurements on azimuthal asymmetries in 
SIDIS 


collaboration 

reaction 

asymmetries 

ref.’s 

HERMES 

e + N -> e + n ± X 

4 sin <t> h a sin 20/, 

UL » A UL 

[66] 


e + N —> e + n ± X 

Asiv, Acoll 

[67] 


e + N -» e + 7T ±f> (K ± )X 

Asiv 

[68] 


e + N -> e + n ± ’°(K ± )X 

Acoll 

[69] 


e + N -> e + n ± (K ± )X 

4 COS 0/, 4 COS 20/, 

*UU » A UU 

[70] 

COMPASS 

p~ 6 LiD -> p-fPX 

Asir , Acoll 

[71, 72] 


p- 6 LiD p-n^K^X A Siv ,A CoU 

[73] 


p-NH 3 -> p-h±X 

Asiv, A Co ll 

[74] 


p~NH 3 —» p^lPX 

Acoll 

[75] 


p~NH 3 — » p~h ± X 

Asiv 

[76] 


p-NH 3 -» p~n ± {K ± '°)X A siv , A Coll 

[77] 


p- 6 LiD -» p-bAX 

4 COS 0/, 4 COS 20/, 

uu » uu 

[78] 

CLAS 

e~ p —> e-jT ±,0 X 

4 sin 20/, 

UL 

[79] 


e~ p —» e~7r°X 

4 sm 0/, 

^LU 

[80] 

JLab Hall A 

e~ 3 He —» e~n ± X 

Asiv, Acoll 

[81] 


e~ 3 He —> e~n ± X 

.cos (0/,—0s ) 

LT 

[82] 


e~ 3 He —> e~n ± X 

.sin(30/,-0 5 ) 

UT 

[83] 


e~ 3 He -> e^X 

Asiv, Acoll 

[84] 


The width has been fitted, the form and flavor dependence 
etc. have been tested. The typical values of the fitted widths 
are e.g. [96], <£~> = 0.25GeV 2 , (P F± ) = 0.20GeV 2 . Roughly 
speaking, this is a quite satisfactory fit. However, it has also 
been pointed out, e.g. in [99] for the TMD FF, that the Gaus¬ 
sian form seems to depend on the flavor and even on z, which 
means that it is only a zeroth order approximation. 

(2) Sivers function: All the data available from HER¬ 
MES [67-69], COMPASS [71-74, 76, 77], and JLab Hall 
A [81, 82, 84] on Sivers asymmetries in SIDIS for pions and 
Kaons have been used for the parameterization. The Sivers 
function is usually parameterized [96, 101-106] in the form 
of the number density f q (x, k ± ) multiplied by an x-dependent 
factor N q (x ) and a k ± -dependent factor li(k ± ), i.e., 

A N fq(x, kj = 2N q (x)h(k ± )f q (x, kj, (5.3) 

where N q (x) is taken as a binomial function of x, 

N q (x. ) = N q x a <{\ - xf*(a q +p q r^/a^ff q \ (5.4) 
and /;(kj_) is taken as a Gaussian, 

h(k ± ) = yf2e(\h\IM x )e- p - IM K (5.5) 

Here the Sivers function A N f q (x, k ± ) is defined via, 

f q /m(x, k ± ) = f q/N (x, k x ) + N f q (x, k ± )!> ■ (p x k ± ), (5.6) 

which is related to the Sivers function fy r (x, k ± ) defined in 
Eq. (3.5) by, 

A N f q (x, k ± ) = -%^x, k ± ). (5.7) 


There exist already different sets such as the Bochum [101— 
103], the Torino [96, 104, 106] and the Vogelsang-Yuan [105] 
fits. One thing seems to be clear that the Sivers function is 
nonzero for proton and it has different signs for u- and d- 
quark, as shown in Fig.2. 



FIG. 2: Example of the parameterizations of the Sivers functions for 
u and d flavors at <2 2 = 2.4(GeV/c) 2 by the Torino group. The figure 
is taken from [104], 

(3) Transversity and Collins function: A simultaneous ex¬ 
traction of them from SIDIS data from HERMES Collabora¬ 
tion [67-70] and COMPASS [71-77] on Collins asymmetries 
in SIDIS and e + e data of Belle [85-87] have been carried out 
by the Torino group [97, 107]. A similar form as that for the 
Sivers function has been taken, e.g.. 

At q(x,kx) = f q (x)[f q (x) + A q(x)] 


Xe- p ^/n(Px), 

(5.8) 

A N D h /q(z,PhT) = 2N q (z)D h i q (z)Kp hT ) 


qi «« 

V 

f-H 

"k, 

1 

X 

(5.9) 


(5.10) 

JVp z ) = A/p z r (1 _ z) <(l_t«2!. 

(5.11) 

Kp hT )- 

(5.12) 


and it has been obtained that also the Collins function is 
nonzero and has different signs e.g. for u — > 7r + or d —* n + , 
as shown in Fig.3. Here, similar to the case for the Sivers 
function, the Collins function A N Di,/ q (z, kf ± ) is defined via, 

Dh/qi(z, Put) = D q/N (z,phT) + ^A N D h / q (z, Phr)s q ■ (k q x p hT ), 

(5.13) 

which is related to the Collins function Hf T (z, pur) defined in 
Eq. (3. 14) by, 

A N D hlq (z, PhT ) = =^H^{z, p hT ). (5.14) 
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FIG. 3: Example of the Torino parameterizations of the transver- 
sity and Collins function. In the left panel, we see the transversities 
xA T q(x) = xhi q (x) for q = u, d\ in the right panel, we see the first 
moments of the favored and disfavored Collins functions. The figure 
is taken from [107], 


The second part concerns the QCD evolution of the TMDs. 
As mentioned earlier, this is a topic that develops very fast 
recently. A partial list of recent dedicated publications is 
[50-63]. QCD evolution equations have been constructed in 
particular for unpolarized TMD PDFs and also for polarized 
TMDs such as the Sivers function. The numerical results ob¬ 
tained from the evolution equations show explicitly that QCD 
evolution is very significant for TMDs. Not only the form of 
the -dependence, but also the width of the Gaussian evolves 
with Q. More precisely, at small k L , Gaussian parameteriza¬ 
tion can be used but the width evolves with Q. At larger k ± , 
the form of k x -dependence is determined mainly by gluon ra¬ 
diation and deviates greatly from a Gaussian and also evolve 
with Q. In Fig. 5, we see an example for the evolution of 
the Gaussian parameterization at small k x \ in Fig. 6, we see 
the evolution of the shape at large k ± . It is also important to 
use the comprehensive TMD evolution rather than a separate 
evolution of the transverse and longitudinal dependences re¬ 
spectively. We show as an example in Fig. 7. 


(4) Boer-Mulders function: It was pointed out that [111] 
the HERMES and COMPASS data on (cos 20) asymme¬ 
try [70, 78] provide the first experimental evidence of the 
Boer-Mulders effect in SIDIS. Studies in this direction has 
been made in [110, 111] to extract Boer-Mulders function 
from the SIDIS data [70, 78] and in [108, 109, 112] to ex¬ 
tract from Drell-Yan data [90-95]. A fit to the first moments 
of Boer-Mulders function of u and d quark is shown in Fig. 4. 
The form was taken again similar to the Sivers function, just 
multiply the Sivers function by a constant, e.g., 

hf(x,k ± ) = A q f^(x,kJ. (5.15) 

However, we would like to point out that the (cos 20) asym¬ 
metry receives twist-4 contributions due to the Cahn ef¬ 
fect [22], A proper treatment of such twist-4 effect involves 
twist-4 TMDs as shown in Eq. (4.39) and in [34]. Because of 
the multiple gluon scattering shown in Fig. 1, the twist-4 ef¬ 
fects could be very much different from that given in [22] the 
results in which corresponds to the case of X. — 1 . A careful 
check might change the conclusion obtained in [108-112]. 




FIG. 4: First extractions of the Boer-Mulders function hf"(x) and 
hf d (x). This figure is taken from Ref.[l 1 1], 

Attempts to parameterize other TMDs such as pretzelocity 
hf T have also been made [113]. Although there is no enough 
data to give high accuracy constraints, the qualitative features 
obtained are also interesting. 



FIG. 5: Example showing the TMD evolution of the Gaussian pa¬ 
rameterization in the low k ± -region. The curves show the evolved 
Bochum Gaussian fits of up quark Sivers function at x = 0.1. This 
figure is taken from Ref. [54]. 

The last thing for TMD parameterizations that we would 
like to mention is the TMD library (TMDlib). We are happy 
to see that, a first version has already been created [114] in the 
year 2014, and updated recently. 


VI. SUMMARY AND OUTLOOK 

In summary, by comparing with what we did in studying 
one dimensional imaging of the nucleon with inclusive DIS, 
we presented a brief overview of our studies on three dimen¬ 
sional imaging of the nucleon with semi-inclusive DIS and 
other semi-inclusive reactions. We summarized in particular 
the general form of the TMDs defined via quark-quark corre¬ 
lators both for TMD PDFs and FFs. We emphasized in partic¬ 
ular on the theoretical framework for semi-inclusive reactions 
at LO pQCD but with leading and higher twist contributions 
consistently. Such theoretical framework is obtained by ap¬ 
plying the collinear expansion technique developed in 1980s 
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i - 0.1, Q - 5 GeV and 91.19 GeV 




M&V) 


FIG. 6: Example showing the evolved k ± dependence in the large k ± 
region. Here we see the up-quark Sivers function at Q = 5 GeV and 
Q = 91.19 compared with the corresponding Gaussian fits at low-k ± 
region at x = 0.1. This figure is taken from Ref. [54]. 



FIG. 7: Example showing the difference between the results of the 
TMD evolution with a DGLAP evolution for x-dependence only for 
unpolarized TMD PDF. This figure is taken from Ref. [55], 


in inclusive DIS to these semi-inclusive processes. We sum¬ 
marized in particular that it applies now also to all processes 
where one hadron is involved. The results obtained in such 
a framework should be used as starting points for studying 
TMDs experimentally. 

At the end, we would like to emphasize that three dimen¬ 
sional imaging of the nucleon is a hot and fast developing 
topic in last years. Many progresses have been made and 
many questions are open. We see in particular that LO pQCD 
leading and higher twists framework for processes where one 
hadron is involved can be constructed using collinear expan¬ 
sions. Factorization theorem for leading twist but with LO and 
higher order pQCD contributions and QCD evolution equa¬ 
tions for unpolarized TMD PDFs and the Sivers functions 
have also been established. Especially in view of the running 
and planned facilities such as the electron-ion colliders, we 
expect even rapid development in next years. 

The overview is far from complete. We apologize for many 
aspects that we did not cover such as the generalized parton 
distributions, the Wigner function, model calculations of 
TMDs, nuclear dependences, and hyperon polarization. 
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